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1. Introduction 

The reconstruction of diagonal elements of density matrix of quantum optical 
states, i.e. of their photon statistics, is fundamental for various applications rang- 
ing from quantum information to the foundations of quantum mechanics j2] 
and quantum optics [Sj. Nevertheless, at the moment, photo-detectors well suited 
for this purpose are not available. Indeed, the choice of a detector with internal 
gain suitable for the measurement is not trivial when the flux of the photons to be 
counted is such that more than one photon is detected in the time-window of the 
measurement, which is set by the detector pulse-response, or by an electronic gate 
on the detector output, or by the duration of the light pulse. One needs a congru- 
ous linearity in the internal current amplification process: each one of the single 
electrons produced by the different photons in the primary step of the detection 
process (either ionization or promotion to a conduction band) must experience the 
same average gain and this gain must have sufficiently low spread. The fulfillment 
of both requisites is necessary for the charge integral of the output current pulse to 
be proportional to the number of detected photons. Few example exist of photo- 
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detectors that can operate as photon counters and each one has some drawback. 
Among these, PhotoMultiplier Tubes (PMT's) |3] and hybrid photo detectors jS] 
present a low quantum efficiency, since the detection starts with the emission of 
an electron from the photocathode. Solid state detectors with internal gain, in 
which the nature of the primary detection process ensures higher efficiency, are 
still under development. Highly efficient thermal photon counters have also been 
used, though their operating conditions are still extreme (cryogenic conditions) to 
be suited for common use [HI Ej . 

The advent of quantum tomography did provide an alternative method to 
measure photon number distributions jSj. However, this scheme, which was ex- 
perimentally applied to several quantum states [S], needs the implementation of 
homodyne detection, which in turn requires the appropriate mode matching of the 
signal with a suitable local oscillator at a beam splitter. Such mode matching, gen- 
erally a not extremely simple scheme to be realized, is a particularly challenging 
task in the case of pulsed optical fields. 

On the other hand, the photodetectors usually employed in quantum optics, 
such as Avalanche PhotoDiodes (APD's) operating in the Geiger mode Ej, 
seem to be by definition useless as photon counters. They are the solid state 
devices, which present the highest quantum efficiency and the greatest stability 
of internal gain. However, they have the obvious drawback that the breakdown 
current is independent of the number of detected photons, which in turn cannot 
be determined. The outcome of these APD's is either "off" (no photons detected) 
or "on" i.e. a "click", indicating the detection of one or more photons. Actually, 
such an outcome can be provided by any photodetector (PMT, hybrid photode- 
tector, cryogenic thermal detector) for which the charge contained in dark pulses 
is definitely below the output current pulses corresponding to the detection of at 
least one photon. Note that for most high-gain PMT's the anodic pulses corre- 
sponding to no photon detected can be easily discriminated by a threshold from 
those corresponding to the detection of one or more photons. 

Nevertheless, recently the possibility of reconstructing photon distribution through 
on/off detection at different efficiencies has been analyzed jll] and its statistical 
reliability investigated in some details ^2]. In addition, the case of few and small 
values of rj ^3] has been addressed. A first experimental application of the method 
of |131 112j has been presented in showing the very interesting potentialities of 
this scheme. In short the method used in ref. is based on the measurement of 
on/off detection frequencies for a certain optical field when varying the quantum 
efficiency of the system, i.e. in practice by interposing calibrated neutral filters 
on the optical path. The beauty of this scheme with respect to alternative ones 
resides indeed in the extreme simplicity that can allow an extensive application of 
it to test optical fields in various applications. 

With the purpose to refine these initial studies, in this paper we present a 
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detailed analysis of the possibility of reconstructing experimentally diagonal el- 
ements of density matrix of quantum optical states and of precisely estimating 
their uncertainties. Various examples will be considered, as PDC heralded pho- 
tons, multi-thermal and attenuated coherent states. 



2. Theoretical methods 

The statistics of the "no-click" and "click" events from an on/off detector, 
assuming no dark counts, is given by 

Po(v) = £(1 - VTQn , (1) 
n 

where g n = (n\g\n) is the photon distribution of the quantum state with density 
matrix g and rj is the quantum efficiency of the detector, i.e. the probability of 
a single photon to be revealed. At first sight the statistics of an on/off detector 
appears to provide quite a scarce piece of information about the state under in- 
vestigation. However, if the statistics about po(r)) is collected for a suitably large 
set of efficiency values then the information is sufficient to reconstruct the whole 
photon distribution g n of the signal, upon a suitable truncation of the Hilbert 
space. 

The procedure consists in measuring a given signal by on/off detection using 
different values r) u (y = 1,...,K) of the quantum efficiency. The information 
provided by experimental data is contained in the collection of frequencies /„ = 
fo(Vv) = nQ U /n u where uq u is the number of "no click" events and n u the total 
number of runs with quantum efficiency rj u . Then we consider expression (^Q) as 
a statistical model for the parameters g n to be solved by maximum-likelihood 
(ML) estimation. Upon defining p v = po(Vu) and A vn = (1 — r] v ) n we rewrite 
expression (JJ) as p u = J2 n A- vn g n - Since the model is linear and the parameters to 
be estimates are positive (LINPOS problem), the solution can be obtained by using 
the Expectation-Maximization algorithm (EM) [T2|. By imposing the restriction 

Qn = 1) we obtain the iterative solution 

K A f 
Qn - Qn 2^ v A rr (*)n { ' 



where f>i/[{f?n }] are the probabilities p u , as calculated by using the reconstructed 
distribution {^«' ) } at the i-th iteration. As a measure of convergence we use the 
total absolute error at the i-th iteration s^> = Y^=o fv ~ Pv[{Qn}] and stop the 
algorithm as soon as goes below a given level. The total error measures the 
distance of the probabilities ^[{gn }], as calculated at the i-th iteration, from the 
actual experimental frequencies. As a measure of accuracy we adopt the fidelity 



= J2 n y g n g^n between the reconstructed distribution and the theoretical one. 



[Author and title] 



4 



3. Confidence intervals 

Our reconstruction of the photon statistics is based on the ML algorithm, which 
by an iterative method reaches the g n best reproducing the experimental data, i.e 
the statistics of 'no-clicks' f v . The confidence interval on the ML determinations 
may be thus evaluated according to the following argument. Due to the fact that 
ML methods are unbiased at convergence we have Pv{q™} = fv If in practice we 
stop at a given iteration L we have fluctuations in the reconstructed {g n } and, in 
turn, Pu{Qn} fv- The fluctuations in the single components g n leads to 

pABu) = {Qn~ Q* n ) (3) 

where we have denoted by g* the true value of the distribution. As a consequence 
we may write bq v n = \p u {Qn) — fu\/A un and, averaging over the different values of 
the quantum efficiency, 

<^ - if L -t (4) 



4. Experimental data 

In the following we describe the reconstruction of photon statistics of several 
optical states. We begin by considering Parametric Down Conversion (PDC) her- 
alded photons. In our set-up, a pair of correlated photons is generated by pumping 
a /3-Barium-Borate (BBO) crystal with a continuous wave (cw) Argon ion laser 
beam (351 nm) in collinear geometry. The crystal has been cut for producing 
type II PDC, i.e. photons with orthogonal polarisation. After having separated 
the photons of the pair by means of a polarizing beam splitter, the detection of 
one of the two by a silicon avalanche photodiode detector (SPCM-AQR-15, Perkin 
Elmer) was used as an indication of the presence of the second photon in the other 
channel, namely a window of 4.9 ns was opened for detection in arm 2 in corre- 
spondence to the detection of a photon in arm 1. The arm 2 "heralded photon" 
was then measured by a silicon avalanche photodiode detector (SPCM-AQR-15, 
Perkin Elmer) preceded by an iris and an interference filter (IF) at 702 nm, 4 nm 
Full Width Half Maximum (FWHM), inserted with the purpose of reducing the 
stray light. More in detail, the output of the first detector started a ramp on a 
Time to Amplitude Converter (TAC) that received the second detector signal as 
stop. The TAC output was then addressed to a multichannel analyser (MCA) 
and a window of 4.9 ns was set by selecting only the respective channels of MCA. 
The quantum efficiency of the arm 2 detection apparatus (including IF and iris) 
was measured to be 20% by using the PDC calibration scheme (see ^E])- Lower 
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Fig. 1: Reconstruction of statistics of quantum optical states. In the highest 
part: Experimental frequencies f u of no-click events compared with reconstructed 
ones (red curve) and theoretical model for the state (black curve). In the lowest 
part, reconstructed probabilities compared with data, a) PDC heralded photon 
(produced in type II PDC cw-regime). b) Strongly attenuated coherent state (He- 
Ne laser beam) 

quantum efficiencies were simulated by inserting calibrated neutral optical filters 
on the optical path. 

A comparison of the observed frequencies f v with the theoretical curve (1 — r] u ) 
is presented in the inset of Fig. 1. In the figure evaluated uncertainties are also 
shown. 

The photon distribution has been reconstructed using K = 34 different values 
of the quantum efficiency from rj v ~ to r\ v ~ 20% with n v = 10 7 runs for 
each r\ v . Results at iteration i = 10 are shown in the figure. As expected the 
PDC heralded photon state largely agrees with a single photon Fock state. The 
estimated relative uncertainty for the p\ element is 0.1 %, confirming the precision 
of the method. One can notice that also a small two photon component and a 
vacuum one are observed. The p2 contribution is expected, by estimating the 
probability that a second photon randomly enters the detection window, to be 
1.85% of px, in agreement within uncertainty with what observed. A non zero po 
is also expected due to background. This quantity can be evaluated to correspond 
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to (2.7 db 0.2)% by measuring the counts when the polarization of the pump beam 
is rotated to avoid generation of parametric fluorescence. Also this estimate is in 
good agreement with the reconstructed po. 

As an example of application of the reconstructed diagonal elements of the 
density matrix, a clear test of non-classical nature of this optical field can then 
be performed by using the evaluated probabilities of finding n photons, p n , to 
estimate the Klyshko [23] parameter K n = (n + l) p ""^" +1 ? which is always larger 
than unity for classical fields. In our case we obtain K\ = (3.2 ± 0.4) • 10~ 4 . 

A further reconstruction of diagonal elements of density matrix has concerned 
a strongly attenuated coherent state. This state has been produced by a He-Ne 
laser beam attenuated to photon-counting regime by insertion of neutral filters. 
Also in this case the same silicon avalanche photodiode detector was used. The 
distribution has been reconstructed with K = 15 different values of the quantum 
efficiency from r\ v ~ to rj u ~ 66% with n„ = 10 7 runs for each r\ v . It agrees 
well with what expected for a coherent state with average number of photons 
|a| 2 ~ 0.02. In the inset of the Fig. 1 the frequencies j v as a function of r\ v are 
compared with the theoretical prediction p u = exp{— 77^ | ct | 2 } ~ 1 — r/ u \a\ 2 . Notice 
that in this case we do not have IF or irises in front of the detector and all the 
other attenuations can be included in the generation of the state: thus the highest 
quantum efficiency is taken to be 66% as declared by the manufacturer data-sheet 
for the photodetector. 

Finally, as a last example, we present here our preliminary data on a single 
branch of PDC in pulsed regime, which is expected to correspond to a multi- 
thermal state 

The state has been generated by pumping a lxlxl mm type I BBO crystal 
by a beam of a Q-switched triplicated (to 355 nm) Neodimium-Yag laser with 
pulses of 5 ns, power up to 100 mJ per pulse and 10 Hz repetition rate. Since 
this source was built with the purpose of having a high spectral selection (as a 
step toward a source useful for realising a quantum memory and quantum logical 
gates based on Kerr effect ;24.) a monochromator with a 0.2 nm selection (ob- 
tained by an entrance 0.025 mm slit) was inserted in the optical path. More in 
details, after having identified the precise direction of PDC emission at 780.2 nm 
by injecting into the crystal a laser beam locked on this wave length at the an- 
gle such to observe stimulated emission, the monochromator was alligned on this 
beam (after a lens) and followed by an objective coupling the transmitted PDC 
to a fiber addressed to photo-detector. The photo-detector was then gated by a 
signal coming from the laser in order to be opened in coincidence with the pulse 
arrival. On/off measurements corresponding to this window were then used for 
the statistics reconstruction, whose results are shown in the figure. 

Also in this case uncertainties in the evaluation of diagonal element of density 
matrix were evaluated according to the discussion of paragraph 3. As shown in Fig. 
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2 the reconstructed statistics of our preliminary data agrees (within uncertainties) 
with the one expected for a multi-thermal case. Nevertheless the fit seems to 
favour a relatively small number of modes (smaller than 10), whilst the number of 
modes dictated by the ratio of the pulse duration (5 ns) and the coherence time 
for the 0.2 nm spectral selection is rather large. If one simply fixes the number of 
modes to 500 or more, the fit becomes worse (the reduced \ 2 increases from 0.5 
to about 10). 

The solution of this problem is postponed to a further data acquisition and to 
a more careful theoretical analysis. One possibility is that the method is reaching 
its limit when the average number of photons is small (smaller than 1 in our case) 
while the number of modes is very large. An alternative solution could derive 
from the presence of a non PDC background, e.g. the introduction into the fit of 
a strong poissonian background can lead to a good fit (reduced x 2 ~ 3 • 10~ 3 ). 

For the sake of completeness, we have also considered the case of multithermal 
distribution with a very low level of photons as well. In this case the optical state 
was a single branch of PDC emission obtained by pumping a type II BBO crystal 
with a cw Argon ion laser beam (351 nm) of 0.3 W. The emission was selected by 
a 4 nm FWHM interference filter (which fixes the coherence time of the field) and 
addressed to an APD single photon detector where on/off counts were collected 
in a 20 ns window. The reconstructed statistics distribution agrees well, within 
uncertainties, with a multithermal state with a large number of modes (see Fig. 
2), but at this low intensity level (0.064 photons on average) the distinction among 
different optical states as coherent, multithermal, thermal or poissonian becomes 
rather weak and no clear distinction emerges from data. 

5. Conclusions 

In summary in this paper we have presented a detailed analysis of the possibility 
of reconstructing experimentally diagonal elements of the density matrix, with an 
accurate evaluation of uncertainties, of quantum optical states by using on/off 
detectors. 

The method that has been used |13|I12| is based on the measurement of on/off 
detection frequencies for a certain optical field when varying the quantum efficiency 
of the system, i.e. in practice by interposing calibrated neutral filters on the optical 
path. 

Various examples have been considered extending and deepening the analysis 
presented in Ref. In particular, we have considered PDC heralded photons, 

multi-thermal and attenuated coherent states. Estimation of uncertainties with 
this technique has been discussed in detail. The results of reconstruction look 
very good for weak coherent and heralded photon state, while a further analysis 
is still necessary for clarifying the case of multithermal state, despite a qualitative 
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Fig. 2: Reconstruction of statistics of quantum optical states. In the highest part: 
Experimental frequencies /„ of no-click events compared with reconstructed ones 
(red curve) and theoretical model for the state (green curve). In the lowest part, 
reconstructed probabilities compared with data, a) multithermal, single branch 
of PDC in nanosecond pulsed regime with 0.2 nm spectral selection. The curves 
correspond to 2 modes, 0.74 average photons, b) multithermal, single branch of 
PDC in cw regime with 4 nm spectral selection. The curves correspond to 10000 
modes 0.064 average photons. 

agreement. 

Altogether our results further demonstrate the potentialities of this technique, 
whose main advantage with respect to alternative schemes resides in the extreme 
simplicity that allows an extensive application for testing optical fields in a wide 
number of applications. 
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